Decoherence by engineered quantum baths 
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We introduce, and determine decoherence for, a wide class of non-trivial quantum spin baths 
which embrace Ising, XY and Heisenberg universality classes coupled to a two-level system. For 
the XY and Ising universality classes we provide an exact expression for the decay of the loss of 
coherence beyond the case of a central spin coupled uniformly to all the spins of the baths which 
has been discussed so far in the literature. In the case of the Heisenberg spin bath we study the 
decoherence by means of the time-dependent density matrix renormalization group. We show how 
these baths can be engineered, by using atoms in optical lattices. 
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Understanding decoherence is central to the descrip- 
tion of the crossover between quantum and classical be- 
haviour Q, Q and it is a crucial issue for the successful 
implementation of quantum information processing 
Although desirable, it is not always possible to fully char- 
acterize the bath and therefore it is necessary to resort to 
ingenious modclizations. Paradigmatic models represent 
the environment as a set of harmonic oscillators Q or 
spin- 1/2 particles In order to grasp all the subtleties 
of the entanglement between the system and its environ- 
ment it would be of great importance to study engineered 
baths (and system-bath interaction) that can be realized 
experimentally and whose properties are amenable of an 
exact solution. Here we discuss a new class of spin baths 
which satisfy these requirements. We show how to realize 
them by means of optical lattices and we find the exact 
solution for the decoherence. We thus provide a way to 
test in the laboratory the emergence of decoherence in 
quantum systems. 

The setup we consider is shown in Fig. ^ Here a two- 
level system (the quantum system) is coupled to a single 
spin of a one-dimensional spin-1/2 chain (the environ- 
ment). The free evolution of the chain is described by 
means of a Hamiltonian which embraces Ising, XY and 
Heisenberg universality classes and it is therefore well 
understood. 

In this context we analyze how decoherence induced 
on the two-level system depends on the internal environ- 
ment dynamics. In particular by varying the parameters 
of the optical lattice we can test what happens when 
the bath enters different phases (critical, ferromagnetic, 
anti- ferromagnetic, etc...). The motivation to analyze 
this model stems from the interest in the decoherence 
due to spin baths both in the absence 0, Q and in the 
presence 0, 0, 0, 0, El of interaction among the 
spins of the bath. In particular Dawson et al. suggested 
that, due to the monogamy of entanglement, decoherence 



depends strongly on the entanglement present within the 
bath 0. The model presented here allows us to frame 
this observation in a broader context, by pointing out dif- 
ferent aspects between the internal entanglement of the 
bath and the induced decoherence of the quantum sys- 
tem. We believe that our analysis may provide a fairly 
general understanding of the relation between decoher- 
ence and ground state properties of spin baths. We will 
first present the exact solution of the theoretical problem 
and then discuss its experimental implementation. 

The Model - The system-bath model shown in Fig. 
is described by the Hamiltonian H = Htl + He + Hin ■ 
The two-level system, characterized by the ground state 
\g) and the excited state \e), has a free Hamiltonian 
Htl = w e |e)(e| and is coupled to the first spin of the 
bath through the interaction Hin = — e|e) (e|crf with e 
being the coupling constant. The Hamiltonian of the en- 
vironment is 
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where af (a = x, y, z) are the Pauli matrices of the i-th 
spin. The constants J, A and A represent the exchange 
coupling, the anisotropy parameter and an external mag- 
netic field. The model defined by Eq. (JIJ has a very 
rich structure 01 . We consider the two separate cases 
of {A = 0, A ^ 0} and {A / 0,A = 0}. When the 
anisotropy parameter is set to zero, and < 7 < 1, the 
model of Eq. belongs to the Ising universality class 
which has a critical point at A c = 1; when 7 = it be- 
longs to the XY universality class. At A = 0, the model 
of Eq. is critical if —1 < A < 1, in the other cases 
the chain has ferromagnetic or anti-ferromagnetic order 
if the anisotropy is positive or negative respectively. Our 
model goes beyond the model of a central spin coupled 
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FIG. 1: a) A sketch of the system plus bath model we con- 
sider in this work. The two level system (at position zero) is 
coupled to the a z component of the first spin of the chain that 
acts as a spin bath. Atoms in an optical lattice can simulate 
this controlled decoherence by means of series of lasers b) and 
displacements of the lattice c) which allow to realize both the 
interaction of the bath with an external magnetic field and 
the anisotropic exchange coupling present in Eq. JQ. 



uniformly with all the spins of the bath and can be more 
easily simulated with cold atoms as discussed towards the 
end of this paper. 

The evolution of the reduced density matrix p of the 
two-level system corresponds to a purely dephasing pro- 
cess. In the basis of the eigenstates \g) , |e), the diagonal 
terms p gg {t) and p e e(t) do not evolve in time. Only the 
off-diagonal terms will decay according to the expression 
p eg (t) = Peg(0)D{t) where D(t) = (\e™t e -i(H T L+H E )t\) 
is evaluated over the ground state of the spin bath. The 
decoherence of the system is fully captured by the func- 
tion £(t) = \D(t)\ 2 sometimes called Loschmidt echo. In 
the following we will study the dependence of L{t) upon 
the spin bath Hamiltonian both in the case of periodic 
and open boundary conditions for a number N of bath 
spins of the order of ~ 10 2 — 10 3 . 

In the case A = and A ^ 0, the function C(t) can 
be calculated exactly. By means of the Jordan- Wigner 
transformation cr+ = ct exp (i7r 5Zfc=i c I c fc); an d a j = 
2cjCj — l, it is possible to map the Hamiltonian of the spin 
bath onto a free fcrmion model which can be expressed in 
the form W E = l^C*, where *+ = (c{ . . . cjy c x . . . cat) 
(cj are the corresponding spinless fermion operators) and 
C = a z <g> A + ia v <g> B is a tridiagonal block matrix with 
Aj.k = -J(Sk,j+i + Sj,k+i) - 2 (A + ej)Sj,k and B jik = 
—J n /(8kj+i — 5j : k+i)- The Loschmidt echo can then be 




0.5 1 1.5 



FIG. 2: a) Loschmidt echo as a function of time for a qubit 
coupled to a N = 300 spin Ising chain; e — 0.25. The various 
curves are for different values of the transverse magnetic field: 
A = 0.5 (black), 0.9 (red), 0.99 (green), 1 (blue), 1.01 (brown), 
1.1 (magenta), 1.5 (cyan). Inset: zoom in at small times, b) 
In the region I, C has a Gaussian decay with a typical scale 
a. Here we plot a/e 2 as a function of A for different values 
of e: 0.01 (circles), 0.025 (squares), 0.05 (diamonds), 0.075 
(triangles up), 0.1 (triangles left), 0.25 (triangles down), 0.5 
(triangles right). The solid line shows the result of a perturba- 
tive calculation at small times, c) In the region II, C oscillates 
around a constant value Coo ■ Here we plot this constant back- 
ground as a function of A for e = 0.25. Points indicate the 
data from the exact solution, while the solid line is the result 
of the perturbation expansion in the coupling constant e. 



evaluated exactly: 

L{t) = det (l - r + re tCt ) (2) 



where r is a matrix whose elements r^j — (\f r ]'\t r j) are the 
two-point correlation functions of the spin chain. Equa- 
tion (J2J) allows us to go beyond the central spin model 
considered so far in the literature and enables us to ad- 
dress explicitly the case of a large number of spins in the 
bath. 

Fig- 12^ shows the generic behaviour of £ as a function 
of time for different values of A (7 = 1). For A < 1 the 
echo oscillates with a frequency proportional to e, while 
for A > 1 the amplitudes of oscillations are drastically 
reduced. As the chain is finite, at long times there are 
revivals, but for N ~ 10 3 there is already a wide interval 
(region II) where the asymptotic behaviour can be an- 
alyzed. It is useful to consider in detail the short- and 
long-time behaviour (regions I and II in Fig. of C. 
At small times the decay is Gaussian, Cj(t) ~ e~ Qt (re- 
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gion I is expanded in the inset of Fig. The scale 

of the Gaussian decay at short times displays a remark- 
able universal behaviour shown in Fig.^b where a(X, e) 
is shown as a function of A. At long times (Region II) 
and for A > 1 the Loschmidt echo approaches an asymp- 
totic value Loo, while for A < 1 it oscillates around a 
value which is constant in time (see Fig. for a qual- 
itative picture). This limiting value strongly de- 
pends on A as shown in Fig. Et. Evidence that Coo de- 
scribes the asymptotic regime can be obtained by com- 
paring the data with the result of an analytical expres- 
sion based on a perturbation expansion in the coupling, 
Coo ~ [l-(e 2 /2) E M0 K^kf l^o)lV(^-^o) 2 ] 2 where 
iV'fc/o) are the excited states (ground state) of He with 
energy E k / . Results obtained from the previous expres- 
sion are plotted in Fig. ^ (solid line) together with the 
exact solution (points). 



i — 1 — 1 — 1 1 1 1 1 1 1 — 1 — 1 — 1 1 1 1 1 1 1 — 1 — 1 — 1 1 1 1 1 1 1 — 1 — 1 — 1 1 1 1 1 




50 100 500 1000 2000 Ji^ 
1 1 I 1 1 1 I 



l(f 



10" 



10' 



10" 



Jt 



FIG. 3: Loschmidt echo as a function of time at the critical 
point A = A c , for different sizes of the chain: N — 50 (black), 
100 (red), 200 (green), 300 (blue), 400 (brown), 500 (violet), 
1000 (orange), 2000 (cyan) and at e = 0.25. The dashed 
line is a fit with a logarithmic decay of the type C c (t) — 
co/(l + ci hit). Inset: Minimum value of as a function 
of N. Numerical data have been fitted with Eq. (dotted- 
dashed line), where l x « 0.997, (3 « 2.369 x 10" 3 . 



Additional information emerge when analyzing the 
scaling of C with the size N of the chain. If the chain 
is not at the critical point and TV is large, both the a 
and the saturation value Coo are almost independent of 
the bath size. At the critical point the situation is rather 
peculiar (In Fig. |3|we plot the decay of C at A c , for dif- 
ferent sizes of the chain) . The decay is very slow in time 
C c (t) ~ In" 1 1, moreover the minimum value C^ reached 
by the Loschmidt echo depends on N as 
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All the properties described so far are typical as long 
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< 7 < 1 (the model belongs to the Ising universality 
class). If the chain is described by the XY model (7 = 0) 
the Loschmidt echo behaves as in Eq. © for A < 1 (the 
model is critical in this range). For A > 1 we found 
C(t) = 1; in this case the coupled qubit does not decohere 
at all. 

In the other case we consider (A ^ 0, A = 0) an ana- 
lytical solution is not available. We then resort to the re- 
cently developed time dependent Density Matrix Rcnor- 
malization Group (t-DMRG) [H 03 with open bound- 
ary conditions. The results shown in Fig. have been 
obtained for chains of N = 100 spins and different val- 
ues of the anisotropy A. We have set a Trotter slicing 
dt = 10~ 2 and a truncated Hilbert space of dimensions 
to = 100. Also in this case the Loschmidt echo is Gaus- 
sian at short times and the scale a of this decay for the 
Heisenberg model is shown in Fig. 0Jd. The dependence 
of a on the different regions of the phase diagram is evi- 
dent. The rate of decoherence is strong in the critical re- 
gion, — 1 < A < 1, because of the presence of low energy 
modes. In the ferromagnetic phase, A > 1, L does not 
decay because the bath is fully polarized. In Fig. 0Jd we 
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FIG. 4: a) Loschmidt echo for a qubit coupled to an N — 100 
spin XXZ Heisenberg chain; e = 0.1 and A = 0.9 (black), 
0.5 (red), (green), -0.5 (blue), -0.9 (brown), -1 (magenta), 
-1.5 (violet), -2 (cyan), -2.5 (orange), b) Decay rate a of the 
Loschmidt echo at small times as a function of A; the dashed 
lines indicate the critical region — 1 < A < 1. 

plot the decay rate a as a function of A. The dependence 
on the critical properties of the chain is again remark- 
able. We notice two qualitatively different behaviours at 
the boundaries of the critical region: at A — +1 there 
is a sharp discontinuity, while at A = — 1 the curve is 
continuous; in the critical region a is constant. 

The rich structure of entanglement in the ground state 
of Hamiltonian suggests the following picture. When 
the decay is quadratic, Jt <C 1, only short-range correla- 
tions in the bath are important and it is therefore natural 
to relate a to the nearest-neighbour concurrence 
We checked this both for the Ising and the XXZ model 
and we found that a is proportional to the concurrence 
for A > 1 for the Ising case and A < — 1 for the XXZ case. 
In the long-time limit long-range correlations become im- 
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portant, in this case the entanglement measured by the 
block entropy 21] or the localizable entanglement [2^ 
seems to be relevant to the decoherence process. 
Experimental implementations - The properties of de- 
coherence discussed so far can be tested experimentally 
in optical lattices. Following the idea of a Universal 
Quantum Simulator described in |23| . the time evolu- 
tion operator associated with H. over a time t can be 
simulated by decomposing the evolution into a prod- 
uct of operators acting on very short times r <C t. 
In our case the operations required are Uj (6) = e l0a ] , 

and Ufu(6) = e l9<T 3 cr ^. For a S {x,y} one can write 



u?? = v; 



Here, Vf = (1 - ia?)/y/2 



are fast homogeneous local unitary operations that can 
be realized with single atoms trapped in an optical lat- 
tice 1^3, each having two relevant electronic levels (\0)j, 
interacting with a resonant laser. These opera- 
tions can be made very fast by simply increasing the 
laser intensity and thereby the Rabi frequency. They 
can be performed either simultaneously on all qubits, by 
shining the laser homogeneously onto all atoms, or se- 
lectively on some of them, by focusing it appropriately 
(see Fig.^i). For our purposes the individual addressing 
is needed only for the atom in position 0, which repre- 
sents the quantum system. Two-qubit operations can be 
performed by displacing the lattice in a state-selective 
way 0, so that state |0)j|l)j+i acquires a phase factor 
e~ ltp ', as experimentally realized in j^. The resulting 
gate Gjj+i((p) can be composed with o~ x rotations to 
yield uj* +1 (0) = e ie [G jd+1 (29)a]a] +1 ] 2 . Since we want 
a different coupling for the {01} pair than for all others, 
we need to erase the effect of the interaction for that spe- 
cific pair using only local operations on atom 0, as in the 
sequences [ofiC/gf(0)] 2 = WW V M? = HUg?(6)} 2 = t. 
To generate each simulation step, we shall first apply 
the — ea z a z term in Eq. (JIJ on all atoms (plus possibly 
other single-qubit terms to simulate the transverse field) 
and then the remainder of two-qubit terms - including 
(e — JA/2)a z a z - while taking care of washing out the 
phases in the (01) pair as described above. An alterna- 
tive scheme |26j , based on tunnel coupling between neigh- 
bouring atoms rather than on lattice displacements, can 
attain the simulation described here for the special case 
7 = 0, but would require some additional stroboscopic 
steps in order to reproduce the general case. Moreover, 
an additional optical tweezer holding atom would be 
needed in this scheme to achieve a system-bath coupling 
different from the intra-bath one. 

Finally we would like to comment on possible exten- 
sions of this approach. It would be quite interesting to 
consider as an engineered bath a three-dimensional op- 
tical lattice. Beside being feasible from an experimental 
point of view, this could be useful in studying for instance 
the situation found in solid-state NMR 271. It would 



be also intriguing to study the Bose-Hubbard model as 
a bath, which would make the experimental realization 
even simpler. Here we chose to focus on spin baths, as 
they are central to the study of decoherence in many 
quantum systems. 
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